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Graph  1 

The  theoretical  foundations  of  normal  mode  theory  and  its  use  are  already 
well  establised.  In  the  future  we  wish  to  consider  range  dependent  problems 
in  terms  of  partitioning,  so  that  the  range  dependent  effects  will  be  corre¬ 
cted  at  the  boundaries  of  the  partitions  and  the  solutions  allowed  to  propag¬ 
ate  further.  In  order  to  apply  this  approach  it  is  necessary  to  find  the 
normal  modes  at  each  boundary.  We  are  interested,  therefore,  in  finding  a 
method  that  will  allow  us  to  calculate  the  normal  modes  in  a  faster,  more 
efficient  manner.  One  of  the  methods  that  holds  promise  of  attaining  this 
objective  is  that  of  transmutation  theory.  Transmutation  theory  is  essen¬ 
tially  a  generalization  of  the  method  of  integral  transformation. 
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Graph  2 

We  are  dealing  at  this  point  with  a  fairly  simplified  model  of  the  ocean 
and  make  the  following  assumptions: 

1)  uniform  depth, 

2)  ocean  is  an  isotropic  medium, 

3)  sound  speed  a  constant 

4)  pressure  p-0  at  the  surface,  and 

5)  ocean  and  its  bottom  are  eventually  (at  some  depth  h) 
underlaid  by  a  rigid  surface. 

Under  these  assumptions  the  excess  pressure  p  satisfies  the  Helmholtz 
Equation. 


prr(r ■ *  ;  pr<r’z)  +  p2z(r'z)  +  k2p<r’z>-  »• 

with  the  boundary  conditions 

p(r ,0)  -  0  and  pz(r,h)  -  0. 

The  second  boundary  condition  (the  boundary  condition  at  h)  is  that  of  a 
rigid  subbottom.  We  solve  in  a  standard  manner  by  using  the  method  of  separa¬ 
tion  of  variables  p(r,z)  -  i/>(z)9(r). 
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THE  IDEALIZED  STURM-LIOUVILLE 

PROBLEM 


THE  IDEALIZED  NORMAL  MODES  SATISFY 


<£"(z)  +  k2  <f>(z)  -  i<f>(z), 

(f)(0)- 0  AND  <f>'( h)  =  0 
WITH  L2  NORMALIZATION  THIS  HAS  SOLUTIONS 


jl.  _  /h  (2m - 1)  7 rz 

rm(z)  “Vo  s,n  5R 


m  =  1,  2,  3, . . . 


i.  -  *  - 


Graph  3 

It  is  well  known  that  for  the  idealized  ocean  with  constant  sound  speed 
the  normal  modes  satisfy  an  idealized  Sturm-Liouville  problem  given  by 


*"<z)  +  k2^(z)  -  i*(z)  , 
and  boundary  conditions 


^(0)  -  0 


*'00  -  0 


It  is  also  well  known  that  this  problem  has  a  complete  set  of  eigenfunctions 
given  by 


<f>  (z) 

m 


(2m  -l)?rz 


m  -  1,  2, 


and  corresponding  discrete  eigenvalues  given  by 


t  -  k 
m 


2  r (2m  -l)i 


m  -  1 ,  2 ,  ... 


These  eigenfunctions  have  been  L2  normalized. 
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We  wish  to  solve  a  similar  problem  associated  with  amore  general  ocean 
model.  The  condition  that  the  sound  speed  be  constant  is  replaced  with  the 
assumption  that  the  sound  speed  is  a  function  of  the  depth  only,  the  other 
assumptions  remain  unchanged.  We  consider  this  model  to  be  the  perturbed 
ocean  model.  The  depth  dependent  equation  gives  a  perturbed  Sturm-Liouville 
problem  for  the  normal  modes  which  satisfy 

0"(z)  +  k2n2(z)V(z)  -  A\Kz) 

with  the  boundary  conditions  expressed  as 

V-(O)  -  0  and  V>'  (h)  -  0. 
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Graph  5 

2 

The  quantity  n  (z)  is  the  index  of  refraction  and  we  represent  it  as 
n2(z)  -  1  +  «s(z) 

where  the  quantity  «s(z)  is  considered  the  perturbation.  We  introduce  the 
parameter  e  which  reflects  the  size  of  the  perturbation  from  its  average 
value.  We  assume  that  the  function  s(z)  is  known  for  the  case  of  interest. 
When  €-0  we  have  recovered  the  idealized  ocean  with  constant  sound  speed. 


5 


TD  7791 


THE  GENERAL  PERTURBATION  APPROACH 


FIND  THE  EIGENFUNCTIONS  AND  THE  EIGENVALUES 
\n  OF  THE  PERTURBED  PROBLEM  IN  TERMS  OF  THE 
EIGENFUNCTIONS  <f>n  AND  EIGENVALUES  Hn  OF  THE 
IDEALIZED  PROBLEM  AND  THE  PERTURBATION  es(z). 


/ 


Graph  6 

In  a  perturbation  approach  we  find  the  eigenfunctions  and  eigenvalues 
Am  of  the  perturbed  (depth  dependent)  in  terms  of  the  eigenfunctions  <f>m  and 
eigenvalues  Jlm  of  an  idealized  problem  (which  we  know)  and  the  perturbation 
€ s ( z )  (which  we  also  know).  Specifically  we  look  for  the  changes  or  correc¬ 
tions  which  must  be  made  to  the  idealized  eigenfunctions  and  eigenvalues. 
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THE  CLASSICAL  APPROACH 


EXPAND  IN  POWER  SERIES  IN  e 


Xm  =  +  eXm  +  e2Xm  +  . . 


AND 

=  <t> m(z)  +  e&)  +  e2  i/$z)  +  . . .  , 

BUT  OBTAIN  EACH  TERM  OF  ^m(z)  ONLY  AS  A 
FOURIER  SERIES 

,(i).  .  0)  ,  .  . 

=  *-  «mp  0p(z>* 

o  =  1 


Graph  7 

The  first  approach  we  examine  is  the  classical  perturbation  approach 
found  in  Titchmarsh.  The  eigenvalues  and  eigenfunctions  of  the  perturbed 
problem  are  expended  in  power  series  of  the  parameter  e . 

The  perturbed  eigenvalue  equals  the  idealized  eigenvalue  plus  corrections, 

\»  -  l*  *  <A<1)  ♦  <2i<2)  +  •••  • 

The  perturbed  eigenfunctions  equals  the  idealized  eigenfunctions  plus  correc¬ 
tions  , 


0m(z) 


-  4>m(Z)  +  €^m(1)(z)  +  «V2)(z)  +  ••  •  • 


However,  we  find  the  corrections  to  the  perturbed  eigenfunctions  only  as  a 
Fourier  series  of  the  idealized  eigenfunctions  and  we  must  find  the  Fourier 
coefficients  a  in  terms  of  the  idealized  eigenfunction,  idealized  eigenvalues 
and  the  perturbation. 
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Graph  8 

In  a  fairly  straightforward  manner  the  explicit  equations  for  the  corrections 
may  be  derived.  We  exhibit  the  corrections  to  the  first  order  in  e.  The  first 
order  correction  to  the  mc  eigenvalue  is  given  by  the  formula 

. (1)  2k2  fh  ,  , r  .  <2m  -l)*zl2  , 

"  TT  J0  s(z>  [  - 2h - J  dz 


The  first  order  correction  to  the  eigenfunctions  is  given  in  terms  of  an  in¬ 
finite  fourier  series  whose  coefficients  are  given  by 


,(1> 

ran 


2k 

h 


2  ,h 


,  ^  .  (2m  -1)ttz  .  (2n  -1)ttz 

s ( z )  sm - ^ - -  sin - ^ - 


dz 


and 


,(1> 

mm 


0. 
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We  say  we  wish  to  transmute  the  solutions  of  an  idealized  Sturm-Liouville 
problem  into  the  solutions  of  the  perturbed  Sturm-Liouville  problem.  A  trans¬ 
mutation  refers  to  a  generalization  of  the  method  of  integral  transformation. 
It  has  the  advantage  that  the  choice  of  the  idealized  system  (the  system  for 
which  the  solutions  are  known)  is  at  our  discretion,  we  have  only  to  find  the 
proper  kernel  function  for  the  integral  transform.  As  a  rule  of  thumb  one 
chooses  the  idealized  problem  to  be  as  similar  as  possible  to  the  perturbed 
problem.  We  choose  for  our  idealized  problem  a  Sturm-Liouville  problem  with 

<T(z)  +  k2^ (z)  -  \<f>( z) 


and  wish  to  transmute  the  solutions  into  solutions  of  the  perturbed  Sturm- 
Liouville  problem  with 

’/'"(z)  +  k2n2  (z)V'(z)  -  AV'(z). 
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Graph  10 


THE  TRANSMUTATION  OPERATION  COMPUTES 

=  <p(z)  +  j  *  K(z,s)  <p(s)  ds 

IN  TERMS  OF  THE  TRANSMUTATION  KERNEL  K(z,s) 
WHICH  MUST  SATISFY 

0  -  0  +  k2  (n*<z)  -  1)  K  =  0. 


We  use  a  transform  of  the  type  given  by 


^  (z)  -  4>  (z)  +  [  K(z ,  s)$ (s)  ds 

m  n  J  h 

where  we  must  specify  the  kernel  function  K(z,s). 

where  the  kernel  function  K(z,s)  must  be  determined.  We  substitute  the  trans¬ 
mutation  into  the  separated  ordinary  differential  equation  for  the  variable 
coefficient  Helmholtz  equation.  This  yields  a  partial  differential  equation 
for  the  kernel  function  K(z,s) 

32K  32K  .  2  .  n 

— r  - r  +  k  £S(Z)K  -  0  . 

dz  ds 
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Graph  11 

We  now  invoke  our  boundary  conditions  for  the  Sturm-Liouville  problem 
with  variable  coefficients.  We  ask  that  ^(0)-0  and  \&'(h)-0.  We  impose  the 
same  bottom  boundary  condition  <£'(h)-0  on  the  solutions  of  the  Sturm-Liouville 
problem  with  constant  coefficients.  Combining  this  with  the  given  transmuta¬ 
tion  we  obtain  two  additional  constraints  on  the  kernel  function  K(z,s) 

~  K(z,h)  -  0  and  2  K(z,z)  +  k^es(z)  -  0  . 

os  oz 

These  two  conditions  and  the  PDE  for  K(z,s)  are  sufficient  to  uniquely  deter¬ 
mine  the  kernel  function 

It  would  seem  that  we  have  made  the  problem  more  difficult,  we  had 
reduced  the  problem  to  an  ODE  and  now  we  must  solve  a  PDE.  Why  not  directly 
solve  the  original  Helmholtz  equation?  The  answer  is  that  we  will  not  dire¬ 
ctly  solve  the  problem  but  instead  will  show  a  way  in  which  the  kernel  func¬ 
tion  K(z,s)  may  be  approximated  as  a  power  series  in  the  parameter  e. 
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"  ' "  '  * 

ONE  METHOD  OF  SOLVING  FOR  THE  KERNEL,  K(z,s): 

LET  M (£,n)  =  K(z,s)  WHERE  2<f  =  z  +  s  -  2h ,  2n  =  z  -  s. 

THEN 

2  Mf„  +  k2(n2(£  +  n  +  h)  - 1)  M  =  0. 

EXPAND  M{£,n)  AS 

oo 

M (£,n)  =  I  k2p  M(pl  (€,n),  WHERE 

p  =  i 

2Mm(<f,/?)  = (n2{a+h)-1)  dcr-J^  W+h)  -  1)d0 

AND 

m<p  +  =  -  Jo  So  M<P’ {a’P)  (n2(a  +  P  +  h>~  dad^- 

- 


Graph  12 

One  method  of  solving  for  the  kernel  function  K(z,s)-  is  to  trans¬ 

form  into  characteristic  coordinates  given  by 

2f  -  z  +  s  -  2h  and  2 >7  -  z  -  s  . 

We  then  get  the  PDE  for  which  is 

2  +k2n2(£  +  n  +  h)M  -  0. 

o 

We  next  do  a  Born  expansion  of  M  in  powers  of  the  parameter  k  .  This  yields  a 
sequence  of  integrals  where  the  first  coefficient  of  the  expansion  is  given  by 

2  (£  ,r?)  -  -f  e  s  (a+h)  da  -[  es(/3+h)  d/3  , 

Jo  J0 

and  succesive  coefficients  in  the  expansion  are  given  by 

M(P+1)(C,r;)  -  -  f  [  M(p)  (a  ,/3)  es(a+/3+h)  da  d/3  . 

J0  J0 
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Graph  13 


Inverting  coordinates  from  (£,r?)  back  into  (z,s)  we  obtain  an  expansion 
of  the  kernel  function  K(z,s).  The  first  term  of  the  the  expansion  is 


f(z+s)/2  f(z-s+2h)/2 

2elC  '(z,s)  -  -  x  i s (a)  da  -kz  «s(a)  da 

Jh  Jh 


and  the  succesive  terms  of  the  expansion  are  given  by 


cK(p+1)(z,s)  -  -k2 


(z-s)/2  p(z+s-2h)/2 


K(p)(a,/3)  ts(a+/3+h) 


We  have  now  obtained  expicit  formulas  for  the  power  series 
kernel  K(z,s)  in  the  parameter  e. 


da  d/3  . 

expansion  of  the 
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TO  SATISFY  THE  BOUNDARY  CONDITION  AT  THE  SURFACE 
<p( s)  IS  CONSTRAINED  BY 

0  *  *(0)  =  <*>(0)  +  |h°  K(0,s)  H s)  ds. 

THEREFORE  WE  DEFINE 

d(z)  =  4>(z)  -  0(0). 

NOTICE  THAT 

d”(z)  +  k2  (d(z)  +  HO))  =  x(0(z)  +  0(0)); 

0(0)  =  0  AND  6  '(h)  =  0. 

ITS  EXPANSION  IS 

d(z)  =  0<°>( z)  +  €0(1)(z)  +  e2e(2\z)  + 


Graph  14 

We  exhibit  the  value  of  the  transmuted  solution  at  the  surface  and  see 
that,  in  general  the  transmutation  does  not  preserve  the  surface  boundary  con¬ 
dition.  Therefore  we  construct  a  new  function 

6{z)-  Hz)  -  HO) 

and  see  that  this  function  satisfies  the  same  boundary  conditions  as  z ) .  We 
then  expand  this  new  function  in  power  series  in  e . 
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Graph  15 

Substituting  6( z)  into  the  idealized  Sturm-Liouville  problem  with  con¬ 
stant  coefficients  and  isolating  the  terms  with  corresponding  powers  of  e  we 
obtain  a  sequence  of  ODE's  each  with  the  same  boundary  conditions  as  the 
original  problem. 

For  the  0  order  terra  it  is  well  known  that  this  homogeneous  ODE  has  a 
set  of  solutions  (normalized  sine  functions)  for  and  corresponding 

eigenvalues  for  A^  ' . 

(0)  We  then  take  the  inner  Products,  0^  with  the  0th  order  equation  and 
6  with  the  first  order  equation  subtracting  we  get  the  expression 

A(1)[  *(0)(s)  ds  -  A(1)-  ^  (A(0)-k2)*(1)(0)  6(°\s)  ds  ; 

J0  Jo 

which  may  be  solved  if  we  know  the  value  of  the  constant  0^(0). 
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Graph  16 

We  can  express  the  perturbed  eigenvalues  to  the  first  order  in  e  as 


A  -  i  +€[i2n_l)xr|  (1)  | 

n  n  L  2h  Nh  n 


and  the  perturbed  eigenfunctions  as 


V2)  -  V2) 


■KM  Ml- 


(2n  - 1 ) 7tz 


'  II  [‘  •  f]-^ 


-l)xz  ,(0)  ] 

2h -  *n  (2)  dz 


[  J^K(1)(z,s)/0)(s)  ds  -  J  [  J  K(1)(z,s)/0)(s)  ds  J  *<°>(z)  dz) 


While  the  correction  to  the  eigenfunction  may  look  complicated  it  is  the  whole 
correction.  We  have  an  explicit  formula  without  having  to  resort  to  Fourier 
series  expansion. 
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rHE  TRANSMUTATION  RESULTS 

(Cont’d) 


\ 


THE  CORRECTIONS  TO  FIRST  ORDER  IN  e  ARE  DEFINED 
IN  TERMS  OF 

4>r)(0)=  Jo  K(1)(0,s)^(s)  ds 

AND 

•  2  /  »-«+2h  \ 

eK'"(z,s)  =  -  ( jh2  |n%|  -  1]df  +  jh  2  [n%)  -  l]dt). 


GRAPH  17 

The  corrections  for  the  eigenvalues  and  eigenfunctions  are  expressed  in 
terms  of  4> ^  '(0)  and  '(z,s).  The  first  constant  is  defined  by  the  expres¬ 
sion 

4>a){ 0)  -  [  K(1)(0,s)/0)(s)  ds, 

n  JQ  n 


and  the  first  order  in  epsilon  term  of  the  power  series  expansion  of  the 
kernel  is  given  by 


€K(1)(z,s)  - 


h 


(I 


(z+s)/2 
[n  (?)  - 
h 


i]  dr 


1 


(z-s)/2  +  h 

[n  (?)  - 
h 


i]  dr 


)■ 


We  again  have  the  perturbation  type  result  where  the  corrections  are  explicit 
functions  of  6,  the  unperturbed  eigenfunction,  i,  the  unperturbed  eigenvalue 
and  the  expression  [n^(z)  -1]  which  is  the  perturbation. 
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CONSISTENCY  OF  THE 
TWO  APPROACHES 

RECALL 

/  —  S 

*  (2n^7T  *  \J\ n%)  -  1]dr  +£2  [„%)  _  1]dfJ  e(0\s)  d8. 

INTEGRATION  BY  PARTS  GIVES 

=  x/0h(n2(f)  ~  [Sin  ---ft  —  df* 


Graph  18 

The  results  of  the  two  methods  look  disimilar,  however  we  can  show  that 
the  results  are  consistent.  To  show  consistency  of  the  first  order  correction 
of  the  eigenvalues  we  start  with 


-P  -  ‘"H  ■ 

We  substitute  for  <^^(0),  including  the  power  series  expression  for  the 
kernel  function.  Integration  by  parts  yields 


which  is  identical  to  Titchmarsh. 
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WE  CAN  SHOW  THAT  THE  FOURIER  COEFFICIENTS  OF 
THE  FIRST  ORDER  CORRECTION  DEFINED  BY 


“nm  -/JVfffc)  «SW  dz- 


ARE  GIVEN  BY 


e“nm  =  h(X^  -  \(n0))  Jo 


-2k2  fh  .  „  ^  .  (2n-1)7rz  ^._(2m-1)7rz 

Si— (5iT  L  (n2W  -  1)  sin  j— g-i— Sin! —  >_  dz 


AND 


eann  ~  °- 


Graph  19 

The  correction  given  by  the  transmutation  is  radically  different  from 
the  classical  results.  However  if  we  express  the  transmutation  correction  in 
Fourier  series  with  respect  to  the  functions  9 '  we  see  that  the  Fourier 
coefficients 


f 

Jn  n 


<1}(z)  /0)(z)  dz, 
n  m 


result.  ...  /i  \ 

After  again  making  substitutions  for  4>K  '  (0)  and  Kv  '(z,s),  several  in¬ 
tegrations  by  parts  and  some  changes  of  variable  these  coefficients  may  be 
used  to  derive  the  classical  Fourier  coeffecients 


h(A(0)-  A<°> 
m  n 


-  f  [n2 

)  J0 


,,  .  (2n  -1)ttz  .  (2m  -l)7rz 

-  1]  sin - sr -  sin - -  dz 


-  0. 

nm 
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